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ABSTRACT 

The MageUanic Stream (MS) is a nearby laboratory for studying the fate of cool gas streams injected 
into a gaseous galactic halo. We investigate properties of the boundary layer between the cool MS 
gas and the hot Milky Way halo with 21 cm Hi observations of a relatively isolated cloud having 
circular projection in the northern MS. Through averaging and modeling techniques, our observations 
obtained with the Robert C. Byrd Green Bank Telescope (GBT), reach unprecedented 3a sensitivity 
of ^ 1 X 10^^ cm~^, while retaining the telescope's 9.1' resolution in the essential radial dimension. 
We find an envelope of diffuse neutral gas with FWHM of 60kms~^, associated in velocity with the 
cloud core having FWHM of 20kms~^, extending to 3.5 times the core radius with a neutral mass 
seven times that of the core. We show that the envelope is too extended to represent a conduction- 
dominated layer between the core and the halo. Its observed properties are better explained by a 
turbulent mixing layer driven by hydrodynamic instabilities. The fortuitous alignment of the NGC 
7469 background source near the cloud center allows us to combine UV absorption and H i emission 
data to determine a core temperature of 8350 ± 350 K. We show that the Hi column density and size 
of the core can be reproduced when a slightly larger cloud is exposed to Galactic and extragalactic 
background ionizing radiation. Cooling in the large diffuse turbulent mixing layer envelope extends the 
cloud lifetime by at least a factor of two relative to a simple hydrodynamic ablation case, suggesting 
that the cloud is likely to reach the Milky Way disk. 

Subject headings: ISMxlouds — ISM:kinematics and dynamics — ISMistructure — methods:data 
analysis — turbulence 



1. INTRODUCTION 

The Magellanic Stream (MS) is our closest and 
most prominent example of a gaseous interaction rem- 
nant. While such circumgalactic structures are pos- 
tulated to represent i mportant sources of fuel for fu- 
ture star formation flKeres et all l2005t iWakker et al.l 
120081: iDekel fc BirnboimI 120061: iBrooks et al.ll2009[ ). the 
mechanisms whereby this material might be accreted 
back into ga.laxies remain unclear. As emphasized by 
iKeres et"all (|2005l ). the multiphase nature of galactic ha- 
los plays an important role during the accretion process 
by modifying gas stripping and infall p rocesses (for ad- 
ditional perspectives on these issues see ISilk et al.l[T987l: 
IGallagher fc SmithI [20051: ITliUmann et al.l 120061) . Being 
close by, the MS offers a unique laboratory to study the 
rate and nature of gas injected from satellites, as well as 
models to assess the fate of the stripped gas. 

The MS trails across much of the southern Galactic 
sky behind the Magellanic Clouds (MCs), passing near 
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the Southern Galactic Pole at about 1/3 of its length. 
It has been shown to have a continuous coherent ve- 
locity tracing across the southern Galactic sky for 140° 
(jNidever et al.ll20To() . just penetrating the Galactic plane 
at some unknown distance from the Galactic center. Di- 
rect distance constraints for the MS are difficult to estab- 
lish much beyond the MCs themselves, which constrain 
the MS head nominally between 52 a nd 61 kpc , based 
on t he MC distan c es es tablished by iKoerweil (|2009( ) 
and iHilditch et al.l ()2005D . With no embedded stars, 
no nearby ionizing sources or stellar absorption sight 
lines yet detected, the MCs remain its only direct dis- 
tance constraint. A geometrical-dynamical analysis by 
IJin fc Lvnden-Belil (2008) found a Galactocentric dis- 
tance that ii icreased to 70 kpc at abo ut 90° along the 
MS length. iStanimirovic et ahl ()2008f ) found a similar 
constraint at this point assuming that the neutral hydro- 
gen (H i) clump size distribution is driven by the Thermal 
Instability (TI). 

Attempts to simulate the large-scale mechanism of for- 
mation and subsequent evolution of this extensive struc- 
ture have focused mainly on tidal interactions with the 
Galaxy and/or ram pressure stripping of the MC gas by 
the Galactic Halo. These simulations had moderate suc- 
cess when it was believed that the MCs have had multiple 
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orbital cycles around the Galaxy (jMoore fc Davisl[T99l 
iMastropietro et al.|[2005l: iConnors et a l. 2006 ~ However, 
global simulations have been made difBcult in recent 
years given the relatively new constraint that the MCs 
are apparently on their first pass through the Galactic 
system (B esla et al . 2007). B esla et al.. (.2010) developed 
a Smoothed Particle Hydrodynamics (SPH) global MS 
simulation in this first-pass scenario maintaining that 
neither tidal or ram pressure mechanisms are as sufficient 
as once thought to stripping the gas. Both ' Besla et al.l 
(jJOlO) and Nidovcr ct al. (2008) propose other mecha- 
nisms to do the heavy lifting, unbinding the gas prior to 
stripping, through tidal resonances between the MCs in 
the former, and in the latter, a gas "blowout" from star 
formation in the LMC region that includes 30 Doradus. 

Whatever the mechanism by which the gas was re- 
moved from the MCs, ram pressure, tidal effects, ini- 
tial turbulent energy, magnetic field structure and sub- 
sequent interactions with the Galactic Halo all may have 
played a role in sculpting the rich filamentary and clumpy 
structure revealed in Figure[TJ This figure shows the scale 
of the entire MS on an H i image produced by the Galac- 
tic All-Sky Survey (GASS, McClure-Griffiths et al. 2009) 
and takes it down in steps to a specific small, but inter- 
esting cloud toward the downstream end which is the 
subject of this paper. The fine scale structure revealed 
in the GASS image an d in oth er rec ent Hi obser v ations 
such as iPutman et aD ()2003aD and iBriins et all (|2005l ) 
is not captured in the N-Body and SPH global models 
mentioned above. Even finer s tructure was revealed by 
iStanimirovic et all (|2002l I2008D in high resolution map- 
ping of the northern tip of the MS using the 305 meter 
dish at Arecibo. These observations revealed extended 
fine filamentary structure and clumps down to the 3.5' 
angular resolution of the telescope. 

Several theoretical studies as well as numerical sim- 
ulations have explored the effect of various hydrody- 
namic instabilities on the formation of small scale struc- 
ture in warm gas such as the MS moving through a 
hot halo (fMori fc Burkert 200l| iQuilis fc M oore 2001; 
iHeitsch fc PutmanI 2009). However, it is still not clear 
how effective these processes are and on what timescales 
they operate. For example, most instabilities appear to 
act on short timescales (~ 100 Myrs) and therefore re- 
quire that the MS is continuously replenished with fresh 
material. The nature and evolution of the small scale 
MS structure (or stripped gas) is also not well under- 
stood. As such material may eventually constitute a sub- 
stantial source for the MW's star formation in the form 
of infalling warm ionized gas (so called "warm drizzle", 
iBland-Hawthorn et al.l l2007f) , observational constraints 
on the physical processes operating on small spatial 
scales in the MS are highly important. 

Numerous observational studies have revea l ed a m ulti- 
phase nature in MS gas. iKalberla fc HaudI ()2006[) and 
IStanimirovic et all (|2008D showed that about 15% of the 
observed Hi clouds have velocity profiles composed of 
warm and cool components, with a velocity FWHM 
of about 25 km and 3-15 km s~^, respectively. 
Matthews et al. (2009) detected the first H i absorption 
lines against radio background sources in the direction 
of the MS close to the MCs indicating spin tempera- 
ture of 70-80 K. Ha measure ments of th e warm ionized 
component in the MS exist (jWeiner fc Williams .1996i : 



IPutman et al.|[2003bD , but with only a few discrete point- 
ings and with the large 1° apertures, providing little 
insight into the processes on the small, arc minute scales. 

While not being able to provide any spatial informa- 
tion, numerous UV and optical absorption studies have 
been crucial for constraining the abund ance of the MS 
gas (jSlavin et al.l 119931: iFox et al.l I20100 . In addition, 
detections o f Ovi absorpti o n from gas associated with 
the MS by iSembach et al.l (|2003D give strong support 
for the existence of an ionized component around the 
MS with T < 10^ K. It is generally interpreted that 
this component represe nts an interface b etween the hot 
Halo gas at T ~ 10^ K ()Fang et al.ll2006D and the cooler 
MS gas. Studies of lower ionization states suggest the 
existence of diffuse envelopes of somewhat cooler, par- 
tially ionized gas that is not visible in the current Hi 
surveys. Specifically, Si ill like ly probes d ifferent phases 
than O VI with T ^'lO^-'*-^ K (jShuU et al.i .2009) and has 
been detected along a sight-line in t he northern MS wit h 
velocities associated with the MS ([Collins et al.l 120091) . 
Recently, iFox et ail ()2010D have performed an analysis 
of many low and high ion species in UV and optical 
absorption, including Sim, against background source 
NGC 7469. From this work, a picture emerges of a dif- 
fuse, multi-phase transition structure between the warm, 
mostly neutral envelope gas detected in Hi and the hot, 
mostly ionized envelope gas dete cted in O vi. 

As pointed out in Stanimirovic etall (poToh . the rich 
multi-phase structure of the MS suggests cloud longevity 
and a slow mass ablation rate. The Kelvin-Helmholtz 
instability (KH) driven by shear flow between cool MS 
gas and the hot Halo is likely to be the dominant mode 
of cool cloud disruption and relatively rapid ablation 
compared to evaporation through pure thermal conduc- 
tion. However, analytical and numerical treatments of 
cool gas in a hot flow suggest that the rapid ablation 
by KH can be moderated by factors affecting the turbu- 
lent mixture in the boundary laye r. These include mag- 
netic fields (Esquivel et al.' '2006 D, re-cooling o f heated 
gas (iBcgclman & Fabian 1990; Kwak et al.l 12011) and 
even ther mal conduction acti ng locally within the mix- 
ture (Vieser fc Hensledl2007al ). 

The focus of this paper is to probe and characterize 
properties of this boundary layer between cool H i clouds 
in the MS and the surrounding Halo gas. This requires 
very deep H i observations, more sensitive than the exist- 
ing Hi surveys. For ex ample, the predomin antly ionized 
component detected bv lCollins et al.l (|2009l ) has the total 
column density of ~ iQi^-is cm"^ and a neutral fraction 
of ^ 0.01. To detect the corresponding neutral gas we 
require H I column density sensitivity of ~ 10^^~^^ cm~^. 
This is at least 5 times lowe r than the most sensi- 
tive H I survey to date; GASS iMcClure- Griffiths et all 
120091) which achieves 3(t column density sensitivity of 
1.6 X 10^* cm^-^. We approach the required sensitivity 
by applying a new method of spatial averaging to char- 
acterize a cloud in the northern MS (shown in Figure [1} 
observed with very deep H i emission spectra obtained 
with the Green Bank TelescopcEl (GBT). Our analysis 

^ The Robert C. Byrd Green Bank Telescope is operated by 
the National Radio Astronomy Observatory, which is a facility of 
the US National Science Foundation operated under cooperative 
agreement by Associated Universities, Inc. 
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Fig. 1. — A journey in scales: From the entire MS, illustrated in the southern Galactic hemisphere map from the Galactic All Sky Survey 
(GASS, McClurc-Griffitlis et al. 2009)), down to the specific cloud studied in this paper. Countercloc kwise from top; High resolution 
GALFA Arecibo Observatory Hi map revealing fine filamentary structure in the northern MS ( Stanimirovic et al.ll2008iV high sensitivity 
GBT map of a section of one of those filaments (this work), and the cloud discovered in that map containing the NGC 7469 sight-line. 



also benefits from the fortuitous location of the NGC 
7469 Seyfert 1 background source near the cloud's cen- 
ter which wa .s extensiv e ly stu died in UV and optical ab- 
sorption by iFox et all ((201C1( ) (hereafter referred to as 
FSWIO). 

Our observing and data reduction strategies are pre- 
sented in Section [2] Section [3] explains the new spa- 
tial averaging method we applied to the cloud to achieve 
< 10^^ cm~^ sensitivity and probe the usually inacces- 
sible cloud periphery and highly ionized cloud envelope. 
We also develop a simple 3D cloud model matched to our 
observations to obtain a robust estimate of cloud's prop- 
erties and gain insight into the internal cloud structure 
and processes acting at its periphery. In Sections 14. H and 
14.21 show that the cloud properties along with results 
from FSWIO strongly suggest an integral WNM/WIM 
core transitioning to a Turbulent Mixing Layer (TML) 
and in Section 14.31 we investigate how consistent the 
cloud's neutral density characteristics are with a photo- 
ionizing environment by comparing the model to the re- 



sults of 3D ionization equilibrium simulations. Finally, 
in Section 14.41 we discuss the lifetime of the cloud and 
assess the possibility of eventually reaching the Galac- 
tic Disk where a portion of its mass could help fuel star 
formation. We summarize our results in Section [S] 

2. OBSERVATIONS AND DATA REDUCTION 

We used the GBT in 2009 and 2010 as part of a larger 
program to map several sample regions in the northern 
portion of the MS. One of these regions is located in the 
northern tip of the MS and was chosen, among other 
reasons, to include the NGC 7469 background source. 
The beam size of the GBT in the 21 cm line is 9.1', 
which corresponds to 18 3 pc at 70kpc, a di s tance con- 
sistent with estima t es by [ jin fc Lvnden-BelH (|2008[ 1 and 
iStanimirovic et all ()2008[ ) which we will assume here- 
after. 

Observations were obtained using the on-the-fly (OTF) 
mapping mode with in-band frequency switching. Con- 
stant Declination (DEC) rows were scanned and stepped 
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every 3.5' in order to oversample the beam by 2.6 in the 
DEC direction. Spectra were dumped at twice that res- 
olution so that beam-spreading in the Right Ascension 
(RA) direction due to scanning remained smah (~ 2%). 
The 36 second target integration time per 3.5' pixel was 
split into four separate scans, which reduced the impact 
of transient signals such as Radio Frequency Interfer- 
ence (RFI) and spectrometer glitches. The GET spec- 
trometer was used with a bandwidth of 12.5 MHz and 
16384 channels, corresponding to a velocity resolution of 
0.161 kms~^ 

The GET data were reduced using a combination of 
the GETIDL and AIPS data reduction packages, and a 
suite of specially designed Interactive Data Language® 
(IDL) programs. To ensure high sensitivity, special 
care was taken during the baseline calibration process. 
First, a modified version of the GETIDL GETFS pro- 
cedure was used to obtain reference spectra to produce 
temperature-calibrated spectra with gross baseline cali- 
bration. Instead of using a frequency-shifted total power 
reference spectrum or a smoothed version of it, as is nor- 
mally done, a noiseless 3rd order baseline model of a 
line-suppressed total power reference spectrum was used 
for reference. This improves the Signal-to-Noise Ra- 
tio (SNR) by without enhancing large-scale baseline 
structure noise which degrades its Gaussian character. 
The line-suppressed spectrum was obtained by selecting 
the minimum of either the reference or signal spectrum 
at each frequency. This "min filter" statistically prefers 
samples without a positive bias from emission of any kind 
and so reduces image artifacts in the final spectra and 
adds a fixed offset to the baseline, easily removed in the 
following steps. The derived reference spectrum was then 
used to perform a {Signal — Reference)/Reference cal- 
ibration of each spectrum. The spectra were then scaled 
by the system temperature and corrected for atmospheric 
opacity using the standard, built-in correction provided 
by GETIDL routines. 

In the next step, the calibrated and roughly baselined 
spectra pass through two additional baseline removal 
phases. Each of these include multiple steps where base- 
lines are modeled and subtracted. In all but one case, 
baselines are modeled as 3rd order polynomial functions 
only after excluding Galactic emission and subtracting 
unblended Gaussian models of any detected emission 
lines in the spectrum. In the first phase, the baseline 
of the scan average for the X polarization is modeled 
and then subtracted from each spectrum in that scan. 
The baselined X emission is then subtracted from the 
Y polarization, effectively suppressing all emission under 
the reasonable assumption in this case that the emission 
is unpolarized. This leaves only the Y baseline, which 
is modeled to 7th order. This special step removes the 
small « 1.5 MHz sinusoidal baseline component present 
in the Y polarization of GET spectra, probably due to 
double transit refiections from the feed to the circum- 
feren tial gap between main rcficctor panels (iFisher et al.l 
[200l . It also removes the average offset introduced by 
the min filter mentioned previously. In the second step, 
each individual spectrum's residual baseline differences 
from the average are removed in two iterations. 

After baselining, spectra from identical pointings were 
combined and interpolated to a regular grid with pixel 



size one fourth that of the sampling interval (or 0.875') 
using the AIPS task SDGRD. This over-sampled the 
telescope beam by a factor of ten in the final data 
cube, providing interpolated resolution to support pre- 
cise averaging along circular paths when spatial av- 
eraging techniques are applied to the cube. Velocity 
spectra were smoothed with a 6th order Hanning win- 
dow, reducing the resolution to 0.966 kms~^. The ve- 
locity range was truncated to —500 to — 200kms~^. 
The rms noise in the final data cube was measured at 
<7t = 3.4 mK by averaging across emission-free veloc- 
ity ranges. This corresponds to column density noise 
<TN = 1.0 X lOi^cm-2 for a 15kms-i FWHM profile. 
This is almost four times better than the 3.8 x 10^''cm~^ 
noise achieved by GASS, the most sensitive survey to 
date (jMcClure-Griffiths et all 120091) when scaled to the 
same FWHM line width of 15kms~"'^. This sensitivity is 
provided at the GET beamwidth of 9.1' compared to the 
Parkes telescope's 14'. 

3. RESULTS: A SIMPLE CORE -I- ENVELOPE CLOUD 
MODEL 

In the data cube we unexpectedly found a roughly cir- 
cularly projected, mostly isolated Hi cloud in the MS, 
located only ~ 2.5' from the direction of the well-studied 
Seyfert 1 background source NGC 7469. Figure [2] shows 
the H I column density (zeroth moment) and velocity 
(first moment, intensity- weighted) images of the cube 
over the velocity range —370 to — 310kms~^. The cloud 
size is ~ 1.6xbeamwidth. After de-convolving the tele- 
scope beam, the intrinsic cloud size is w 1.2 x beamwidth. 
The cloud is therefore just resolved. The Hi column 
density at the cloud center is w 5 x 10^^ cm~^ and 
although it trails off to the north, east and south, it 
connects to the cloud complex to the west. The veloc- 
ity image shows that this connection involves an abrupt 
~ —10 km velocity shift from the ~ —340 km s^^ 
gas at the cloud's apparent western edge to the neigh- 
boring ~ —350km s^^ component labeled "A" of the 
western complex, which is mixed in with a nearby ~ 
—345 km s~^ component labeled "E". Note also that the 
cloud itself has a distinct ~ -f-5km s~^ gradient from 
west to east. The abrupt positive velocity shift from the 
western complex along with the continued positive gra- 
dient across the cloud suggests it may have broken off 
from the complex and is becoming entrained in the sur- 
rounding ambient gas. 

Figure [3] shows the Hi spectrum at the center of 
the cloud. The emission clearly has a component at 
w -340kms-i with a FWHM w 20kms-^ However, 
an even lower level, broad component lies at slightly 
more negative LOS velocities as made evident by the two- 
component Gaussian fit overlaid on the spectrum. The 
width and weakness of the wider envelope line causes 
us to consider the possibility that it is an artifact of the 
baseline removal process, but we reject that as there is in- 
dependent evidence of components associated with both 
the core and envelope line velocities in the UV absorption 
spectra along the nearby NGC 7469 sight line reported 
in FSWIO (discussed in Section [3|). 

To enhance sensitivity of our observations and char- 
acterize how properties of the two components evolve 
outward from the cloud center, we introduce the spatial 
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Fig. 2. — The Hi column density (a) and velocity (b) images of the region containing the NGC 7469 sight line and the Hi cloud that it 
passes through. The velocity range is —360 to — 320kms~^ and the velocity image column density threshold is 2 X 10^* cm~^. The images 
show that, although the large cloud complex to the west has a similar velocity range, the cloud is connected to it only by diffuse bridging 
gas and is well isolated in velocity. The cloud complex region was excluded from the circular path averaging as indicated. 
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Fig. 3. — Hi spectrum at the cloud center reveals two components, a ^ 20kms ^ wide component at —340 km s ^ on top of a weaker, 
s 60kms~^ wide component. Individual Gaussian fits to each component are shown (dashed lines) as well as the composite fit (solid line). 



averaging method and apply it to the data cube in Sec- 
tion [3TTJ This produces a profile of average column spec- 
tra versus projected angular distance, p from the center 
with increasing sensitivity, well past its apparent periph- 
ery. In Section 13.21 we purge the average profile of com- 
ponents that are unassociated with the cloud. Finally, in 
Section [331 we develop a simple parametric 3D spherical 
model of the cloud optimally matched to the purged 2D 
projected profile to obtain a robust estimate of the prop- 
erties and structure of the cloud versus physical distance, 
r from its physical center outward. 

3.1. Spatial Averaging Method: Raw cloud Profile 



To perform spatial averaging of H i spectra we start 
with the following assumptions: (i) the cloud is roughly 
circularly symmetric, (ii) the cloud symmetry extends 
past its apparent boundaries, and (iii) along each line 
of sight on a circular annulus from the cloud center, H i 
gas has similar properties. The small velocity gradient 
of ~ -1-5 km s^^ on the eastern side of the clump (Fig- 
ure[2Ib), discussed above) shows that the symmetry as- 
sumption is compromised, but the gradient is relatively 
small with respect to the width of both cloud compo- 
nents identified in the central spectrum of Figure [31 Its 
effect on the results of the averaging process that follows 
are considered in Sections 13.21 and T 
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We start frora the pixel with the highest H i brightness 
teraperature (see Figure [2]) and average H i spectra spa- 
tially along circular annuli of increasing projected radius 
p. As a result, we can plot the averaged Hi spectra as 
a function of projected angular distance from the cloud 
center as shown in Figure EJa). The most important 
advantage of this spatial averaging method is that it re- 
duces the noise without compromising angular resolution 
in the essential radial dimension, which contains most of 
the spatial information in the nominally circularly pro- 
jected cloud. The noise improvement in the averaged 
spectra increases with p due to an increasing averaging 
path length, thus compensating for decreasing emission 
from the cloud center. Assuming infinite resolution and 
uncorrelated, zero mean noise at each pixel, the measure- 
ment error for a circular cloud theoretically decreases as 
1/ ^2'Kp. In reality, error improvement will be limited 
to the extent that the noise is zero mean and uncorre- 
lated from pixel to pixel. The normal gridding process 
increases pixel-to-pixel correlation, while residual struc- 
ture from imperfect baseline modeling violates both of 
these criteria and will generally be the limiting factor in 
most applications of this method. While we apply the 
method here to a roughly circular cloud, we note that 
the same method can be applied on filamentary struc- 
tures by assuming linear symmetry and averaging along 
lines parallel to the filament main axis. 

The result of the spatial averaging method is shown 
in Figure lllja). On the x-axis we plot the averaged Hi 
spectra at a given projected distance p from the cloud 
center, shown on the y-axis. The first spectrum at p = ' 
corresponds to the single H i spectrum through the cloud 
center (shown in Figure 3), and is dominated by the main 
emission component centered at w ~ — 340kms~^. As p 
increases, the primary emission component weakens and 
the broader, underlying component strengthens as the 
sensitivity improves. There is also a distinct blue- ward 
velocity gradient in the broad component spectrum as p 
increases. The distinct narrow and broad emission com- 
ponents of the cloud are suggestive of a spherical struc- 
ture with a (kinematically) warm core and a warmer en- 
velope: with increasing p, the line of sight passes through 
less of the core and more of the envelope and eventually, 
the envelope dominates through limb brightening. We 
explore this further through modeling in Section 13.31 

In addition to emission at the main component veloc- 
ities, there appear to be one or more weak lines in the 
spectrum that are well separated from the main emis- 
sion {e.g., at —250km s""'^) and therefore assumed to be 
unassociated with it. The evolving profile also contains 
components at similar velocities, but unrelated to the 
cloud's central components as the widening circle of av- 
eraging passes through clumps of emission that are spa- 
tially separated on the sky. These are also assumed to 
be unassociated. In the next step, these unassociated 
components are purged from the projected profile. 

3.2. Spatial Averaging Method: Profile Line Modeling 

Using a semi-manual process aided by a graphically 
interactive program written in IDL, all Hi spectra in 
the averaged profile were modeled as multiple, blended 
Gaussian functions. The process is as follows: Each H i 
spectrum of the profile is viewed in order of increasing 



p and initial estimates of all potential Gaussian compo- 
nents are provided by the user to a blended Gaussian 
fitting routine. Each component of the resulting fit is 
then checked and any with SNR < 2 are removed with 
the remaining components used as initial conditions for 
a second fit. At this point, since Gaussian decomposition 
is not unique, the user views the fit and can accept it or 
try again with new initial estimates. 

The result of this modeling process is shown in Fig- 
ure m^b). While many unrelated lines pop in and out 
with increasing p, there are two main Gaussian compo- 
nents that maintain continuity from p = ' outward and 
so are considered to be the components of the cloud. All 
others are to be subtracted. It is important to stress 
that the purpose of this line modeling step is not to ac- 
curately model all the emission in the profile but only to 
identify and subtract components in the observed profile 
that are likely to be unrelated to the cloud. The resul- 
tant "cleaned" profile, which contains only the two ob- 
served primary components directly associated with the 
cloud, is shown in|4l^c). Note that there is a prominent 
discontinuity in the cleaned profile which occurs at the 
point where the broad component becomes undetectable 
as part of the blended fit underneath two rising, but unre- 
lated components that dominate it. This is an artifact of 
imperfect line modeling. Where the broad component is 
too weak to contribute its own line model to the blended 
fit, the other components incorporate it into their line 
models and then it gets subtracted. 

Another limitation of the line modeling is cross- 
coupling, where a stronger component affects the param- 
eters of a weaker overlapping component. These effects 
are not easily quantified and could dominate any un- 
certainties due to noise alone. Imperfections notwith- 
standing, the lines associated with the two components 
of the cloud are modeled as part of the cleaning pro- 
cess and provide a rough model of the cleaned profile 
as shown in Figure SJd). Figure [S] shows how the mod- 
eled line width, central velocity, and H i column den- 
sity parameters of these two components evolve with p. 
. Figure [Sfa) shows that the two components have dis- 
tinctly different line widths. The narrow core component 
has a relatively constant FWHM ~ 20 km s~^ while the 
broader envelope has an apparently large outward gra- 
dient from ~ 50km s^^ at cloud center to ~ 80km s~^ 
at p ~ 20' where the line falls below threshold. Al- 
though the gradient may be exaggerated by parameter 
cross-coupling with the core component at low p and with 
subtracted components at higher p, it could be significant 
and is physically consistent with turbulent and/or ther- 
mal broadening as it approaches the hot ambient Halo 
medium. 

The central velocity profiles of Figure [5{b) illustrate 
the parameter coupling concerns mentioned previously. 
At small p the stronger core component may "pull" the 
velocity of the weaker envelope, while the roles are re- 
versed at p ~ 9 ', where the core weakens and its velocity 
drops sharply toward the envelope's. The sharp drop in 
core velocity and the step change in envelope velocity are 
probably not entirely real. The core velocity is likely con- 
stant for the most part at ~ —336 km with the sharp 
drop clearly a modeling artifact since this is not possible 
given the beam smoothing. The envelope's apparent neg- 
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Fig. 4. — Spatially averaged spectral profiles in the line modeling phase of the analysis where the measured profile is "cleaned" of emission 
not related to the cloud: (a) Observed spectral profile of the cloud showing the intensity spectrum (horizontal axis) averaged along circular 
paths of increasing distance from the center (vertical axis), (b) Blended Gaussian emission line model of the entire profile, (c) Cleaned 
spectral profile of the cloud where Gaussian emission line models of all except the lines associated with the two components of the cloud's 
central region are subtracted is used to constrain the optimally matched 3D parametric cloud model in the next step, (d) Lines associated 
with the two cloud components. For image display purposes only, profiles have been smoothed from 0.966 kms~^ to 4.83kms~^. Note 
that the vertical striations in the noise are the consequence of a noise correlation length of roughly four pixels due to interpolation. 
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Fig. 5. — Cloud core and envelope Gaussian component line model parameter profiles, (a) Line Width (FWHM); (b) Velocity; and (c) 
Column Density. Error bars reflect parameter errors due to noise only. Other effects such as cross-coupling between components (see text) 
are not accounted for. 



ative velocity gradient may be exaggerated by "pulling" 
from the core at low p but is large enough that it could 
be at least partially real, approaching ~ —350 km s^^ at 
the periphery. 

Gradient or not, there is a difference in the central 
velocity between the two components of between 4 and 
14km with the core component clearly lagging the 
envelope component along the LOS. This is very inter- 
esting because the large negative velocity of this cloud 
indicates infall and a cometary morphology with the core 
leading the envelope is expected. Although the morphol- 
ogy may be hidden in projection or in the noise of Fig- 
ure [H the core in the spatially integrated profile is clearly 
lagging in velocity. This is also apparent in Figure [3] 



(the central cloud spectrum) indicating tha t this is not 
an art ifact of the spatial averaging process. 'Briin s et al.l 
(|200H ) similarly modeled the line profiles of HVC125-h41- 
207 along simple perpendicular cross sections (i.e., no 
spatial averaging) and demonstrated that its projected 
on-sky head-tail morphology had a cold core leading a 
warmer envelope. Similarly, many other HVCs with a 
head-tail structure have a lagging tail, which is usually 
interpreted as belonging to gas stripped and slowed down 
through interactions with the surrounding medium. This 
interesting, apparently anomalous velocity structure of 
our cloud is addressed in Section SI Some of the gradi- 
ent in line width and central velocity of Figure [5][a) and 
(b) could be caused by the velocity gradient asymmetry 
toward the eastern side of the clump. We discuss this in 
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Section [ 

Figure [Sjc) shows the well-behaved Hi column density 
of the two components. The core component is roughly 
Gaussian with FWHM ~ 10', only slightly larger than 
the 9.1' beamwidth. This reflects the barely resolved 
core. The envelope component has an almost constant 
column density of ~ 2 x 10^^ cm~^. It is worth emphasiz- 
ing that the broad envelope component which has a sig- 
nificant H I column density even at large p, constitutes a 
large portion of the cloud's total neutral mass and prob- 
ably an even larger portion of its total mass since it is 
likely more highly ionized than the core. In Section 13.31 
we quantify this in the context of the 3D cloud model 
developed there. 

As a by-product of the component fitting process the 
residual noise of each averaged H i spectrum is measured 
after subtracting the line fit. This is plotted in Figure |6] 
as a function of projected angular distance p from the 
cloud center. This figure clearly demonstrates the ex- 
tremely high sensitivity of this spatial averaging method, 
where cjn reaches below 4 x 10^^ cm~^ at about 10 ' from 
the cloud center. Also plotted is the theoretical noise as 
discussed in Section[321 which also incorporates the noise 
correlation effects of interpolation. Measured noise does 
not improve as rapidly as it could, probably because of 
limited resolution and residual baseline artifacts at low 
p and the correlation introduced by the gridding pro- 
cess as mentioned before. Note that at the beamwidth, 
p ~ 9', correlation effects are minimal and the noise im- 
provement approaches the theoretical limit. However, at 
p > 12', the averaging path begins to be truncated by 
the excluded region and then by the borders of the cube 
and so the noise flattens and begins to increase. 

3.3. 3D Cloud Model Fitting 

We now use the two-dimensional projected profile of 
Figure W^c) to investigate a simple three-dimensional, 
core-envelope cloud model supporting our observational 
claims regarding the cloud structure. Since the cloud 
is not spatially resolved well in our measurements, sub- 
tle spatial structure in the model is not appropriate, 
but since spectral resolution is very high, subtle spec- 
tral structure can be incorporated. We define a spheri- 
cally symmetric cloud model with three spatial and one 
velocity dimension as: 



nHi{r,v) ■- 



2-K<T^{r) 



, Tc < r < rc + re 
,r>rc + re 



(1) 



where 



(Tc(r) = (jyco + m„cr, 
<Je{r) = (Jveo + m^er, 



r = \/a;^ + 2/^ + 



(2) 
(3) 
(4) 
(5) 



and cloud spatial coordinates x, ?/, z are defined so that 
the z-axis is along the LOS, x is east to west on the sky, 



while y is north to south. The cloud consists of a core 
with a constant H i volume density Uco and has a radius 
Tc- Beyond the core is the cloud envelope with a constant 
density rieo which extends up to a radius of r^- At the 
cloud center, r = 0, the core and envelope have a LOS 
velocity Vco and v^o, respectively. The spectra of the core 
and envelope are Gaussian functions centered at Vco and 
VeO, respectively. The velocity dispersion cr„ of the core 
is allowed to linearly vary with r, but the central velocity 
is held constant since the observed column profile (Fig- 
ure [5|) shows little variation. Significant gradient in the 
velocity and dispersion of the envelope is suggested by 
the profile, so both are allowed to vary linearly with r. 

Note that this model incorporates a radially dependent 
velocity scalar field that is restricted to the LOS compo- 
nent of the vector field. As such, the model, although 
very sensitive to the directly measured LOS motions as- 
sociated with a cometary feature seen in projection, is 
insensitive to radial motion from the cloud's 3D center. 
So, expansion or contraction effects cannot be inferred 
from the model. Parametric variables used in this model 
are defined in Table [TJ 

In practice, the function of Equation ([T|) was used to 
generate a model grid with three spatial dimensions plus 
one for LOS velocity. The spatial grid pixel spacing was 
the same as that of the observed radial profile. The grid 
encompassed the entire spherical cloud which spans twice 
the profile's radial extent in each dimension, but was ex- 
tended in the x and y dimensions by two beamwidths 
at zero density to allow for beam smoothing. The spa- 
tial averaging method was applied to the 3D model by 
summing the spectra in the z direction at each (x, y) 
pixel to obtain a spectral profile to be compared with 
the observed cloud profile. The weighting used was the 
inverse of the measured profile variance at each pixel. 
Prior to comparing simulated and observed profiles, the 
model array was processed to simulate observing effects. 
This involved approximate beam smoothing with a 2D 
Gaussian having FWHM equal to the GET beamwidth, 
truncated at ±FWHM. An equivalent averaging of spec- 
tra along circular paths was then applied as was done in 
the case of the observed data cube. 

We used the AMOEBA algorithm ()Nelder fc MeadI 
|1965[) to find the model parameters that minimize the 
weighted least mean-squared error of the 3D model's sim- 
ulated observed profile. AMOEBA finds the local mini- 
mum in the parameter space that is nearest to the initial 
guess. All of the parameters in this model relating to 
the velocity spectrum shape are easily estimated from 
and well constrained by the observed profile data. They 
are also relatively independent of the core and envelope 
sizes. For these parameters, the error surface should be 
well behaved with a single minimum. The core and enve- 
lope sizes and densities are highly inter dependent, and 
the error surface would have a very broad minimum given 
that the beam smoothing is comparable to the profile's 
scale width. Also, there may very well be multiple min- 
ima. Because of these concerns, a sweep of the initial 
guesses for these two size parameters was performed in 
order to find the lowest local optimum over a range of 
values. For each combination of initial core and envelope 
size guesses, approximately self-consistent initial density 
values were derived from the observed column density 
at p = and p ^ Tc, respectively. The optimization 
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Fig. 6. — The measured RMS noise of the Hi column density for a 15km FWHM Hne (solid line) versus the projected distance 
from the cloud center. Also shown is the theoretical error improvement (dashed line) assuming zero mean, uncorrelated noise and infinite 
resolution. Measured noise does not improve as quickly at first partially because the effects of finite resolution are more important for short 
averaging paths. Note that close to p ~ 11 ' measured noise approaches the theoretical value of 2.4 X 10^'' cm~^, and then increases due to 
path truncation from exclusion of the western region of the cube. This is a very deep measurement and maintains error in the periphery 
below 4.0 X 10^^ cm-2. 



TABLE 1 

Core plus Envelope Model Opimized Parameters 



Parameter 


Description^ 




Value 




Core normalized radius (arcmin) 




7.32 ±0.37 


ricO 


Core number density at center ( '^^^'^ cm 




(5.27 ±0.41) X 10-3 


fcO 


Core velocity (km s"-"^) 




-336.3 ±0.2 




Core velocity dispersion at center (km s~ 




8.80 ±2.94 


rria-c 


Core velocity dispersion slope (ISS-^- — ) 




0.20 ±0.61 


Te 


Envelope normalized thickness (arcmin) 




18.51 ±0.41 


neO 


Envelope density ( ^^j^^'^cm^'^) 




(8.83 ±0.26) X 10"" 


VeO 


Envelope initial velocity (km s"-"^) 




-349.6 ±0.4 




Envelope velocity slope ('^"^ °. — ) 

^ ./ i V arcmin ' 




-0.441 ±0.264 




Envelope initial velocity dispersion (km s 




25.8 ±3.0 


rrio-e 


Envelope velocity dispersion slope ( °. — ) 

^ .7 r- r- \ arcmin ^ 


0.157 ±0.253 



^ Scalable units use D, the LOS distance in pc. 



was run at these combinations and the solution with the 
minimum error at optimum was chosen. Optimized pa- 
rameter values are shown in Table [T] along with la er- 
ror estimates. The parameter error vector Ap = (Ap„) 
was estimated by taking the square root of the diagonal 
elements of the covariance matrix [J'^WJ]"!. Here, 
W = [l/c,^,!] is the diagonal weighting matrix used in 
the optimization, am are measured error vector elements 
(flattened matrix of rms noise residuals from line mod- 
eling), and J = [Jmn] — [^^] is the Jacobian matrix at 

the optimum solution (determined numerically). 

Figure [7] shows the cloud modeling results. Figure[7{a) 
is the spatially averaged observed projected spectral pro- 
file purged of all but the main emission component (same 
as in Figure EJc)), used as the optimization target for 
the AMOEBA algorithm. Figure [7Jb) is the projected 
spectral profile of the optimized 3D cloud model without 
beam smoothing effects. Figure ^c) is the same as (b) 
with beam smoothing, which is used to compare with 
i.e., subtract from, (a) in the optimization process. Fig- 
ure [Tfd) shows the residual from this subtraction. The 
residual structure at p from ~ 20 to 24 ' is a consequence 



of the optimization "smoothing" out the discontinuity 
in the cleaned envelope emission at the points where it 
meets line detection threshold. The low residual levels in 
the image show that our simple core-envelope 3D cloud 
model, when similarly projected onto a 2D spectral pro- 
file, clearly compares well with the observed data. This 
supports our interpretation of narrow and broad veloc- 
ity components tracing internal structure of a single Hi 
cloud. 

The quality of the fit for the model is quantified by 
comparing the RMS of this modeling residual, amodei — 
6.17 mK (calculated between —450 and —250km s~^, 
the range of emission) , to the RMS noise of the observed 
profile, CTnoise = 6.20 mK measured in the previous line 
modeling step. This represents a Vk^ode; " ^no»sel = 
0.7 mK RMS modeling error, only 11% of the noise. 
The simplicity of this model and the low residual error 
strongly suggest that this is a reasonable representation 
of the observed data. We therefore proceed with further 
data analysis under the core-envelope cloud model. 

The model parameters are summarized in Figure [51 
where values of H i volume density, central line veloc- 
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(□) "Cleaned" Spectral Profile 
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Fig. 7. — Cloud modeling results: (a) Cleaned spatially averaged spectral profile purged of non-cloud components, (b) The projected 
spectral profile of the two-component 3D cloud model without beam smoothing, (c) The projected model profile with beam smoothing, 
(d) The model profile residual error, or the difference between the cleaned profile and the beam-smoothed model profile. For image display 
purposes only, profiles have been smoothed from 0.966 km s~^ to 4.83km s~^. 



ity and line width for the two components are plotted 
as a function of radial distance, r from the cloud cen- 
ter. The best-fit values for the Hi volume density of the 
cloud core and envelope are 5.3 x 10"'^ and 8.8 x 10~* 
cm~^, respectively, scaled to a reference distance of 70 
kpc. For compar i son, using an independent method 
iStanimirovic et all (|2002[ ) estimated the H i volume den- 
sity for several MS clouds in this part of the MS at 
~ 5 cm~^ assuming a distance of 60 kpc. The density 
of the modeled core scaled to this distance translates to 
6.1 X 10~^ cm"'^, which is ~ i/s of that estimate, but the 
former is based on clumps with > 5 times the central col- 
umn density so a higher neutral fraction is to be expected 
due to more effective self-shielding. The LOS velocity of 
the model core is constant, as constrained by the model, 
and the nominal velocity of the envelope gradually de- 
creases from —344 to —351km s~^. The error bars in- 
dicate considerable uncertainty in this gradient which is 
only ~ 10% of the line width. 

As mentioned previously, the small observed asymme- 
try in the observed velocity field in Figure Elb) could 
contribute to the velocity centroid and width gradients 
observed in Figure [S] (and modeled in Figure [S]). We 
can estimate the magnitude of this effect by quantifying 
the observed velocity gradient asymmetry in the core, 
assume that it continues into the envelope, and then es- 
timate its affect on the model's averaged gradients. From 
Figure [2l[b) , ~ +5 km s^^ gradient occurs over approx- 
imately one beamwidth, or « 10' resulting in a velocity 
gradient of ~ +0.5 km s~^arcmin^^. Since the asym- 
metry is present for less than one-half of the averaging 



path, the combined effect is about one-half the measured 
gradient, or « -1-0.25 km arcmin"^. Referring to Ta- 
ble [11 note that this could account for most of the ob- 
served core line width gradient to^c which is dominated 
by the uncertainty due to noise anyway, so the effect of 
the asymmetry on the core line width gradient is insignifi- 
cant. The model did not include a core velocity gradient 
because there appeared to be little variation, but this 
positive gradient asymmetry may have compensated for 
a roughly equal and opposite gradient amounting to only 
w 2 km s~^ across the « 7' radius of the core. We also 
note that the negative envelope line width gradient ma-e 
could be partially compensated by the asymmetry, but 
the « -1-0.25 km s~^ arcmin""'^ effect is comparable to the 
error due to noise and so not of significant concern. 

Table [2] summarizes the global properties of the model 
cloud in more physically relevant terms useful for discus- 
sion in Section [H Those properties that are distance- 
dependent are scaled for a reference distance of 70 kpc 
LOS. The units indicate how the quantities scale at other 
distances. It is particularly noteworthy that while the 
envelope has about a factor of two lower H i column den- 
sity, it is quite extended and is seven times more massive 
than the core. 

4. PHYSICAL INTERPRETATION OF THE MODEL 

Using the highly sensitive analytical technique based 
on a simple 3D parametric model described in Section [3l 
we have obtained observational evidence of the cloud 
having a distinct neutral core with a WNM-like veloc- 
ity line width, surrounded by an envelope with a much 
wider line width, described by the model parameters of 



Transition from Magellanic Stream to Halo 



11 



-335 : 
-340 : 
-345 - 
-350 - 
-355 - 



- Core 



-Ifflllll 



- Core 
" Envelope 



0.2 0.3 0.4 0.5 
- (kpc • D/70 kpc) 



(a) Component density 



0.0 0.1 0.2 0.5 0.4 0.5 
r (kpc • D/70 kpc) 

(b) Component velocity 



r (kpc • D/70 kpc) 

(c) Component line width 



Fig. 8. — Characteristics of the optimized two-component core-plus-envelope spherically symmetric cloud model: (a) Stepped uniform 
density profile as imposed by the model, (b) The velocity profile shows an abrupt blue-ward shift at the core-envelope transition followed 
by a moderate blue-ward gradient in the envelope, suggesting an outward motion of the core with respect to the envelope, but both are 
only ~ 15% of the envelope line width and therefore, not very significant, (c) Line width profile with the core consistent with a WNM 
phase, and the envelope ~ 9x kinematically warmer. The line width gradients are not significant given the errors, but consistent with 
warming and/or increasing turbulence toward the periphery. These small velocity and line width gradients may also include the effects of 
a small asymmetry in the velocity structure of the cloud (see text). 



TABLE 2 
Model Cloud Properties 



Parameter 


Description'' 




Value 


[l,b] 


Galactic coordinates (deg) 




[83, -45] 




LOS to Galactocentric angle, D > 50 kpc 


(deg) 


< 9.4 




Core radius (j^^pc) 




147 


fie 


Mean Core neutral density (^^^^EHcm"^) 




5.3 X 10-3 




Mean Core velocity dispersion (km s~^) 




10 


Vc 


Mean Core LOS velocity (km s"-*^)) 




-336 


NhIc 


Central Core Hi column density (cm~^) 




4.8 X 10^® 


MhIc 


Core Hi mass {{^^)2Mq) 




2.5 X 10^ 


Te 


Envelope thickness ( YO^pc 




380 


fie 


Mean Envelope neutral density (^^2^P£cm 




8.9 X 10"* 


^ve 


Mean Envelope velocity dispersion (km s~ 




29 


fie 


Mean Envelope LOS velocity (km s~^)) 




-358 


Nnie 


Central Envelope Hi column density (cm" 




2.1 X 10^* 


Mnie 


Envelope Hi mass ((Tof^)^^©) 




1.8 X 10* 


Mhi 


Total Hi mass ((,^)2Mq) 




2.1 X 10* 


MHmin 


Minimum total H mass'' ({t^^)'^Mq) 




4.1 X 10^ 



Scalable units use D, t he LOS distance in pc. 
b Assumes Hii/Hi > 19 (|Fox et al.ll2"010l ') 



Figure [51 We now delve into physical interpretations of 
the observed spatial and kinematic structure. 

The most interesting and puzzling aspect of our re- 
sults so far is the negative velocity offset of the envelope 
with respect to the core which ranges from —5 km 
at the center and gets larger further out to —10 km s~^ 
implying that in this infalling structure, the envelope, 
presumably consisting of gas components stripped from 
the core and entrained by the relatively static Halo, is 
falling faster, leading the core. We would expect a typ- 
ical cometary feature to develop with entrained, enve- 
lope gas trailin g behind the core with a lagging veloc- 
ity as found bv lBriins et al.l (|2001[ ) in a structural anal- 
ysis of HVC125-I-41-207 previously mentioned in Sec- 
tion 13.21 Their analysis of highly resolved observations 
of the HVC revealed WNM-like gas in a cometary enve- 
lope surrounding a well-resolved CNM-like core. In con- 
trast, our cloud's barely resolved WNM core corresponds 
to their entire cloud and our envelope corresponds to 
possible boundary gas surrounding their HVC, well be- 



low their detection threshold. Still, we might expect an 
analogous structure in the boundary gas, in transition 
to the Halo. Our cloud shows no clear projected on- 
sky cometary feature, but given the high LOS velocity it 
would probably be hidden in full projection behind the 
cloud. However, this orientation would directly show the 
lagging envelope velocity signature of a cometary feature 
yet we see it clearly leading. 

The most likely explanation for the leading envelope is 
suggested by the velocity structure of Figure HJb) which, 
as discussed in Section [3l indicates the cloud may be a 
large clump of entrained gas originating from the west- 
ern complex and now lagging it. The detected envelope 
velocity range of ~ —350 km s~^ to ~ —345 km s~^ cor- 
responds closely to the velocity range of the two nearby 
components (A and B) of the western complex. The de- 
tected envelope therefore, may be an extended, diffuse 
structure that both the cloud and the western complex 
are embedded in. Although of low Hi column density, 
^ 2 X 10^^ cm^^, it is quite extensive and may contain 
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a significant amount of H i. 

However, our spatial averaging method provides fur- 
tfier constraints on tlie nature of the diffuse envelope. 
Since the region excluded in the spatial averaging in- 
cludes most of the extended structure and bridging gas 
to the west past ~ 14' from the cloud center, our 3D 
cloud model is not influenced by this western region. We 
can therefore assume that although the projected circular 
symmetry assumption of the averaging method is weaker, 
the core-envelope cloud model (Section 3.3) reasonably 
represents properties of the envelope gas in the vicinity of 
the cloud, suggesting a physical association between the 
cloud core and its immediate envelope. Therefore, we 
proceed with interpretation under the assumption that 
the core is embedded in the surrounding envelope, but 
that the envelope is part of a larger structure in the local 
hierarchy closely associated with the core. 

In addition to our Hi observations, the UV spec- 
troscopy of NGC 7469 by FSWIO provides rich comple- 
mentary information about the cloud. Numerous absorp- 
tion lines have been detected from both low and high 
ionization species at MS velocities. FSWIO concluded 
that neither photoionization nor single-temperature col- 
lisional ionization can explain the observed column den- 
sities of high ions. This suggested the existence of highly 
multi-phase plasma, with a cooler region traced by Si iv 
(peaking at 10^'^ K), and a hotter region being traced 
by Ovi (peaking at 10^-^ K). On the other hand, the 
detected low ionization species of O i, C ii and Si ii trace 
WNM and WIM gas with a temperature < 10** K as does 
Hi. 

Table [3] summarizes the parameters of our observation- 
based H I model components (calculated along the same 
sight line) and measured components of three of the low 
ions from FSWIO. The core component central velocities 
line up quite well within a few km s~^ of the H i and the 
wider component central velocities are also quite close, 
but have a wider range. The velocity dispersions, ay 
of the core components are all comparable, but those of 
the envelope velocity are widely scattered. We interpret 
these data in the following sections. 

4.1. Cloud Envelope 
4.1.1. Smooth conduction- dominated boundary layer 

The detected envelope Hi line has an average mod- 
eled velocity dispersion aye = 29 km s~^ (Table[5]) which 
does not vary significantly across its full extent, which is 
2.6 Tc (Figure[5]). Although the observed velocity FWHM 
might suggest a slight increase away from the cloud cen- 
ter, this is not well constrained with our angular resolu- 
tion. This dispersion is ~ 3 times higher than what is 
expected for the WNM/ WIM gas at - 10* K and could 
be interpreted as being due to gas at kinetic temperature 
T < 1.0 X 10^ K with an additional turbulent component. 
One possible explanation is that the envelope represents 
a conduction-dom inated (evaporat i ve) bo undary layer, 
along the lines of iCowie fc McKei ()1977D . in relatively 
smooth transition from the ~ lO^'K WNM/ WIM core to 
the ~ lO^K WIM Halo. We do not c onsider radiative 
effects on the interface as analyzed by iMcKee fc Cowid 
([1977) since the estimated core radius of 147 pc (Ta- 
ble [2]) is well below the critical radiation radius, Rrad 
defined in that analysis as the radius above which radi- 



ation effects become significant. The analysis, based on 
cooling at solar metallicity yields Rrad > 1.6 kpc for an 
assumed final (Halo) temperature T > 10^ K and density 
n < 10-''cm-2 (jSembach et al.ll2003[ ). This is already 
much larger than our estimated cloud radius. The lower 
metallicity in this part of the MS of ~ 0.1 solar (|Fox et al.l 
120101 ) implies a lower cooling rate. Thus the mean free 
path in the boundary layer, represented by the scale 
length A(T), increases (jMcKee fc Cowid ITQTTI Equation 
7) making for an even higher critical radius. This makes 
the non-radiative assumption even more comfortable. 

Interpreting the envelope as single-phase gas in a 
conduction-dominated boundary layer, we first estimate 
the maximum detectable gas temperature that would 
leave sufficient neutral fraction to produce the observed 
2.1 X 10^^ cm^^ neutral envelope column density (Ta- 
ble [2]) through a realistic column depth. A rough equi- 
librium calculation was performed for Hydroge n using 
an a pproximate collisional ionization formula (jPraind 
120111 Equation 1 3.11) and int erpolated recombination 
rate coefficients (jSpitzeii 119981 Table 5.2) under opti- 
cally thin conditions. It showed that envelope gas of 
density < 0.011 cm^"^ (twice the core neutral density of 
Table [2]) at a temperature > 3.8 x 10'* K, would be so 
highly ionized as to require a column > 140 kpc (twice 
the assumed distance) to obtain the measured column 
density. So, the observed envelope would have to trace 
gas no warmer than 3.8 x 10* K , with kinetic (thermal) 
dispersion ak < 18 km s^* . Given the kinetic compo- 
nent, the turbulent component required to produce the 
total dispersion would have at > 23 km s^^. 

Although no conduction-dominated analyses or simu- 
lations that we are aware of consider intrinsically turbu- 
lent cloud structure, they predict the temperature pro- 
files of static (evaporativ e/condensing) boundary layers. 
iDalton fc Balbuj (119930 considered both classical and 
saturated heat flux and derived temperature proflles as a 
fu nction of r/vc and the global saturation parameter, ag. 
In lGnat et all (|20Tol) Fi gure 1, these profiles are plotted 
for a range of ao values. In our case, we assume Thim = 
Tuaio > 10^ K and for our estimated maximum envelope 
temperature, we have T/Thim < 3.8 x 10*/10^ = 0.038 
where r/vc gets only slightly larger than one, regardless 
of saturation parameter. This would make the transi- 
tion to our maximum envelope temperature on order of 
0.01 Tc. In fact, even if the above detectability analysis 
is ignored and assuming the envelope dispersion is com- 
pletely thermal, with corresponding T^ax = 1-0 x 10^ K, 
the transition distance is still only on order of 0.1 Tc. 
This is clearly inconsistent with our extensive detected 
envelope which reaches a distance of 2.6 Tc beyond the 



iVieser fc Henslerl (|2007bD performed simulations that 
included classical and saturated heat flux (their model 
R3) for a cloud of similar total mass to ours and under 
ambient temperature and density on the same order as 
the Halo. Although nearly sonic flows and a thin turbu- 
lent layer (~ 5% of Tc) developed in these simulations, 
no systemic turbulence developed nor was it considered. 
The initial cloud radius, at 41 pc, is less than 1/3 the size 
of our cloud, but wc assume that things roughly scale 
with radi us and ambient temperat ure as do the derived 
profiles of IDalton fc Balbui (|1993l ). A temperature pro- 
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TABLE 3 

Low Ion Line Properties toward NGC 7469 



Core Envelope 



Line (eV) (km s"!) (km s^^) (km s^i) (km s^^) 



Hi 21 cm'' 13.6 

Oi A1302.168'^ 13.6 

Sin A1260.422= 16.3 

Cii A1334.532'^ 24.4 



-336 ±0.2 9.7 ±5.0 

-337 ±1 5.6 ±1.4 

-335 ±1 7.1 ±1.4 

-330 ±1 9.2 ±1.4 



-357 ±0.4 28 ±5 

-364 ±1 4.2 ±4.2 

-368 ±2 9.2 ±2.8 

-360 ±1 16.3 ±2.8 



^ Ionization potential to next level. 

^ This work, mean model values 2.5 ' off-axis (toward NGC 7469) with 
9.1 ' beam. 

FSWIO Table 1, b values scaled to cr„, aperture is 0.2 x 0.06 ". 



file is not provided, but if we interpret the density profile 
of their Figure 12 assuming roughly constant pressure, 
the steep drop in density from r/vc ~ 1 to about 10 
times that of the ambient Halo at r/vc ~ 1.2 suggests 
the temperature at that point would be T ^ 10^ K for 
Tfiaio ^ 10^ K. From this temperature profile, our de- 
tected envelope gas should extend no further than 20% 
of core radius (again assuming that the velocity disper- 
sion is purely the r mal). This is much further than the 
iDalton fc BalbusI ()1993[ ) profiles predict, but still less 
than an order of magnitude as deep as our observed 
envelope, which extends as far as 260% of core radius, 
even after convolution with t he G ET beam is taken into 
account. IVieser fc HensleJ (j2007b) also simulated the 
cloud with self-gravity (their model R4) and then again 
with cooling/heating added (their model R5), with each 
compressing the transition range further. We can there- 
fore state that our observed envelope extends further, 
by more than an order of magnitude, than simulated or 
derived results predict for an evaporating cloud with- 
out systemic turbulence. The impact that turbulence 
throughout our core and envelope would have on these 
predictions is not known, yet one would expect it to be 
significant. However, an order of magnitude increase in 
the extent of the low temperature range of the boundary 
layer seems unlikely. 

4.1.2. Turbulent mixing layer 

As an alternative interpretation, the large velocity 
dispersion of the cloud envelope may be the result of 
the superposition of multiple, more or less independent 
clumps of WNM/WIM gas turbulently intermixed in 
the boundary layer with warmer phase gas along the 
li nes of the Turbulent Mixing Layer (TML) described 
in lBegelman fc FabianI (|1990D . This approach provides a 
simpler and more supportable explanation for our obser- 
vations. In this picture the broad velocity width of the 
envelope is due to the superposition of discrete clumps 
of ~ 10"* K WNM gas (as found in the cloud core. Ta- 
ble [2]), photoionized by the background and Galactic 
UV radiation, and turbulently dispersed in a mixture 
with much warmer gas approaching HIM temperatures 
with T - 10^ K. The superposition of muhiple WNM 
('^ 10'' K) clumps, each with a^c ^ 10 km s~^ turbu- 
lently dispersed with at ^ 27 km s~ will produce an 
approximate Gaussian shape with cr„e ^ 29 km s^^, sim- 
ilar to the observed total envelope dispersion. 



The low metal ion envelope velocity dispersions of Ta- 
ble |3] should also be consistent with this picture. Oi and 
Sin envelope components have comparable dispersion, 
but are very low compared to Hl Cii has significantly 
higher dispersion than Oi and Sin, yet less than Hl 
It is important to emphasize that Hi is observed with 
a 9.1' beam and averaged all around the cloud's vicin- 
ity through the spatial averaging we applied, so it would 
sample a large number of these WNM clumps over a wide 
range of dispersed velocities. On the other hand, the 
metal ions are observed with a "pencil" beam through 
the center of the cloud, passing through a much smaller 
number of the clumps and having a smaller dispersion. 
The larger C n dispersion could be explained by its sig- 
nificantly higher ionization potential allowing it to sur- 
vive into warmer and more turbulent gas surrounding the 
WNM clumps. 

The TML model, first suggested by 
iBegelman &: FabianI ()1990f ). explains turbulent mix- 
ing of cool and hot gas arising from the KH or shear 
instability induced by the velocity difference across 
the interface between the two types of gas. The TML 
consists of a continuum of phases between the cooler 
10"* K (WNM/WIM) stripped gas and intermediate 
temperature gas at T = ^/TcooiTuot ~ 10^ K (assuming 
Thot = Tnaio ~ 10^ K) coexisting in a turbulent mixture. 
Also, since cooling is very efficient at this intermediate 
temperature, they suggest that some of the ablated gas 
would rejoin the cooler phase rather than be lost to the 
hot phase g possible means for slowing the ablation 
process and extending cloud lifetime. 

Hydrodynamic simulations of WNM gas clouds passing 
at supersonic velocities through the HIM provide further 
qualitative support to our observations. In simulations, 
WNM clouds consistently develop a hierarchy of smaller 
clumps peeling off the cloud as a result of KH insta- 
bility, si milar to the picture o f WNM/WIM clumps we 
propose (lEsouivel etall |200l IVieser fc Hensleil l2007al: 
iHeitsch fc PutmanI l2009t iKwak et all 120111 ). In the 
present case, our cloud core is probably a very large 
clump recently separated from the western complex of 
Figure [2] as discussed above. Another consistent fea- 
ture of the simulations is the development of vorticular 
flow along the shearing interface of the main cloud as 
well as clumps and protuberances that develop there. A 
most striking exam ple are the early evolution simulations 
shown in Figure 2 of lVieser fc Hensleil ()2007aD where vor- 
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tices are fully developed in the non-conductive case and 
appear to be forming in the more slowly developing con- 
ductive case . In th e longer-term multi-ion simulations of 
iKwak et alJ (j2011t ) the material is tracked as the develop- 
ing TML entrains it. The simulation represented in their 
Figure 2 extend to much larger times and the cloud has 
properties more like our cloud than the s i mulatio ns of 
lEsauivel et all (|2006[ ) or iVieser fc Hensleii (|2007aD . As 
the gas falls behind the HVC it mixes with the ambi- 
ent gas producing intermediate temperature gas. The 
temperature and velocity structure of the ablated ma- 
terial increases in complexity and range with time. At 
later times, the simulations indicate a wide 10'^ to 10^'^K 
temperature range in the ablated material mixture with 
'-^ 40 km of velocity range in the ~ 10^ K gas near 
the cloud axis. This is similar to the observed velocity 
dispersion of the H i envelope. Als o significant , although 
not easy to discern in Figure 2 of iKwak et al.l ()2011[ ). is 
the vorticular transport of gas from the mixing layer to 
the central region where it reverses direction with respect 
to the entrained gas as discussed by the authors. They 
also describe it as a means by which gas cooled in the 
mixing layer replenishes the cooler gas at the center as 
suggested by iBegelman fc FabianI ()199ClD and mentioned 
above. This vorticular motion in the simulations is con- 
sistent with our proposed picture of clumps dispersed in 
both leading and lagging velocities with respect to the 
cloud in our observed envelope. 

4.2. Cloud Core 

The observed H i core component, with its mean line 
width consistent with T < 1.2 x 10"* K, is reasonably 
interpreted as the neutral component of a central con- 
centration of moderately turbulent WNM to WIM sup- 
plying cooler gas to the TML after being stripped from 
its periphery. The kinetic temperature T and turbulent 
dispersion at of the gas can be estimated from the line 
velocity dispersions, cri and a2 of two co-spatial species 
with atomic masses mi and m2 by solving the two equa- 
tions: 



hT 

mi 

hT 

m2 



^2 J2 

a, = an. 



(6) 
(7) 



More than one species could be incorporated by finding 
an optimum solution for a set of equations for all co- 
spatial species, of course. 

The components of the three metal species at or near 
the H I core velocity listed in Table [3] should sample a 
line of sight through the WNM core and the (photoion- 
ized) WIM envelope that probably surrounds it. How- 
ever, they may not all be strictly co-spatial. The density 
profile of each species can be quite diflFerent along the line 
of sight due to different ionization potentials responding 
to a decreasing UV fiux toward the shielded core center. 
H I and O i have the lowest ionizing potentials and so 
survive mostly toward the center and concentrate there. 
Si II and C ii have much higher potentials and, surviv- 
ing the higher fiux in the periphery, are distributed more 
evenly throughout the core. Unless gas along the line 
of sight has a uniform temperature and turbulent dis- 
persion out to a distance where Si ii and C ii are mostly 



depleted, these two ions would be biased with respect 
to Hi by sampling gas near the likely warmer and more 
turbulent periphery near or even into the TML. H i and 
O I have nearly identical ionization potential which leads 
to nearly identical ion fraction profiles, sampling points 
similarly along the line of sight, heavily weighted toward 
the self-shielded center of the cloud. As a result, we can 
use Hi and Oi to separate the thermal and turbulent 
components and be confident that the result refiects the 
mean along the line of sight, weighted heavily toward 
the center of the core and not affected significantly by 
the TML. 

Using the velocity dispersion values for H i and O i 
from Table [3] and solving the above equations yields T — 
8350 ± 350K and turbulent dispersion at — 5.3km s~^, 
both of which are quite consistent with expectations for 
WNM/WIM. It is interesting to note that if we assume 
that Si II and C ii are at this same temperature through- 
out the core column, their turbulent dispersion compo- 
nents come out to 6.9km and 8.9km s~^, respec- 
tively. This is consistent with Sin, with its ionization 
potential higher than H i, surviving further into an in- 
creasingly turbulent periphery and C ii, with an even 
higher potential, surviving further still into even higher 
turbulence as the WNM/WIM core transitions to the 
TML. 

4.3. Cloud Ionization 

Our picture of the cloud is of a partially ionized sta- 
ble WNM/WIM core, surrounded by a TML, in which 
diffuse WNM/WIM clumps are dispersed. These small, 
dispersed clumps are exposed to not only ionizing extra- 
galactic and Galactic radiation, but are also embedded 
in the ionizing radiation produced locally by the coUi- 
sional processes in the warmer phases of the TML. We 
therefore assume that the clumps are highly ionized and 
essentially transparent to Galactic UV and extragalactic 
background, which reaches the core at full strength. We 
now investigate whether a 3D photoionization model of 
the cloud core exposed to the Galactic and extragalac- 
tic UV radiation field can reproduce a neutral compo- 
nent consistent with our core model. We use a version 
of the 3D Hyd r ogen- only ionization equilibrium code of 
iWood fc Loebl ()2000D . The code was set up to include 
a Galactic and extragalactic ionizing fiux. A 3D spher- 
ical cloud of uniform volume density is discretized onto 
a linear Cartesian grid with 129 cells on a side. The gas 
is assumed to be isothermal, with recombination rates 
based on a kinetic temperature of 8000 K. 

The gas can be made clumpy with a fractal-generated 
density m odulation calculated using the algorithm of 
Elmeg recn (1997) ( as described in se veral papers e.g., 
Mathis et al. 2002; Wood et al.l [20051) . This algorithm 



leaves a fraction of the gas as the smooth component and 
redistributes the remaining gas into hierarchical clumps. 
The fractal dimension / and parameters Ni through 
govern the process. For more on these parameters and 
their effect on density structure, see the above-mentioned 
references. In this study, our purpose is only to include 
the general effect of dumpiness on the density profile 
shape for no other reason than to add some measure 
of realism to the cloud. It is not to emulate large-scale 
structure (the original purpose of this algorithm) nor any 
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particular ISM environment. We empirically adjusted 
the f smooth and fractal parameters of the algorithm to 
produce, after ionization, only one distinct neutral struc- 
ture representing the cloud core (with some parameter 
settings, there were two cores) while still producing dis- 
tinct substructure. We settled on f smooth — 0.85, fractal 
dimension f — 2.6 with A^i through set to 64, 16, 16, 
16 and 16, respectively. 

In order to demonstrate that photoionization alone can 
produce a neutral core consistent with our observations, 
we exposed a spherical cloud of uniform average density 
to ionizing radiation appropriate to our assumed distance 
of 70kpc. A density of 5.5 x 10~^ cm~^ was assumed 
which is somewhat higher than that of our core model 
(Table[2]), thus allowing for a moderate ion fraction at the 
center. Again this is merely a demonstration illustrating 
what the actual profile, unresolved by our observations, 
might look like based on exposure to a realistic ioniz- 
ing field. The ionization fiux used is a combination of 
constant isotropic extra ealactic ionizing backgr ound of 
7.9 X lO^cm-^s-i from iFaucher-Gigufere eFall (j2009j!| 
and directio nal Galacti c flux of 1.3 x 10^cm~^s~^ from 
the model of iFox et a l.l (j 2005r p| calculated at the cloud's 
Galactic coordinates from Table [2] and assumed LOS dis- 
tance of 70 kpc. The results are shown in Figure O To 
obtain the ionization results, we kept the density and 
fractal parameters constant and adjusted the cloud ra- 
dius in order to obtain the observed central H i column 
density of Nc = 4.8 x 10^* cm"^ (Table [2]) as measured 
along the Galactocentric LOS, defined in the local rec- 
tilinear coordinate system as the z-direction. When cal- 
culating average density profile and column density, the 
cloud center is determined from the 3D centroid of the 
neutral density distribution. 

Figure shows the resultant neutral density struc- 
ture in the x — z plane after the cloud has reached ion- 
ization equilibrium. Figure Wijo) shows the cloud H i and 
total Hydrogen density profile through the centroid along 
the LOS. The photoionized cloud component is more 
extended than the neutral core. Note the asymmetry 
with the neutral core to right of center. This is due to 
the directional Galactic radiation component, providing 
increased flux along the LOS, which is aligned closely 
with Galactic Center. The remaining neutral core has 
a FWHM size of 319 pc. Figure IHl^c) shows the average 
radial Hi and total Hydrogen density profiles from the 
simulation. Also shown is our observation-based model 
core H i profile for comparison (from Table [2]) . 

The neutral density profile (Figure injb)) of the ionized 
cloud core approximates a composite exponential profile 
dropping rapidly near the more neutral center and then 
a longer scale takes over as it becomes more ionized at 
larger radii. This is typical of ionization fronts, as in 
iZheng et al.l (|2002| ). where logarithmic plots show three 
slope regimes, a flat range if the core is fully shielded, 
followed by a steep region with moderate ionization and 
then a shallower one with high ionization. Since in our 
example the core is not fully shielded, there is no flat 
region. Our crude core model of constant density de- 
rived from observations approximates this more realistic 

Calculated from energy spectrum provided online by the au- 
thors at www.cfa.harvard.edu/~cgiguere/UVB.html 

^ From data kindly provided by Joss Bland-Hawthorn. 



density profile as a constant, roughly equal to the cen- 
tral density truncated at about half the full extent of the 
core, which has a radius of 319 pc. 

This phoionization model also provides information 
about the neutral fraction of the cloud. We calculate the 
neutral fraction from Figure |9jb) by comparing neutral 
and ionized hydrogen profiles and find a fraction of 24% 
along the central LOS. This is three times higher than 7% 
determined by FSWIO based on CLOUDY simulations 
for low ions. FSWIO also determined a lower limit for the 
incident ionization parameter, U — n^/nn > 5 x 10~*. 
The equivalent U for our example is9.7xl0~^, a factor 
of 5 lower. To get the FSWIO value of [/, our simulations 
must be at an either unrealistically low or high distance 
(< 20 kpc or > 200 kpc) to get high enough Galactic 
or low enough njj, respectively. Realistically, the ioniz- 
ing fiux and therefore U would be even lower if partial 
shielding of the extragalactic fiux by the Galactic disk 
and the neighboring western complex were taken into ac- 
count, increasing the discrepancy. Given the insight we 
now have into the structure along the NGC 7469 sight 
line, it can be argued that for the core and envelope 
structure we propose, FSWlO's CLOUDY analysis over- 
estimates both U and the ion fraction for the low ion core 
environment. Their analysis used the total metal column 
densities, which include significant metal ion column den- 
sity from envelope velocity components, as targets for 
the simulation. If these absorption components occur in 
highly ionized WIM fragments dispersed in the TML as 
we propose, their contributions would not be well repre- 
sented as part of an integral CLOUDY slab model. Also, 
higher potential ions like Si ii and C ii would dispropor- 
tionately sample warmer portions of the ambient colli- 
sional mixture. Furthermore, local TML collisional FUV 
radiation could contribute to the ionization state of the 
WIM fragments, further biasing the modeling results. It 
would be interesting to repeat the FSWIO simulations 
using only the metal column densities associated with 
the core velocity as targets. Given the high slopes in 
FSWIO Figure 5 (left) for the high potential ions, this 
could result in a significantly lower U and ion fraction 
more in line with our results. 

In summary, we showed that we can reproduce the ob- 
served H I column density and the rough size of the cloud 
core relatively easily and with simple assumptions, in the 
scenario where the core represents the remnant of a larger 
cloud exposed to Galactic and extragalactic background 
ionizing radiation. The phoionization exercise also sug- 
gests that a large percentage (76% in this example) of 
the warm phase core of this cloud is ionized. 

4.4. Cloud Lifetime 

If we assume that the mixing layer is merely a transport 
mechanism for the WNM/WIM of the cloud to ablate 
into the Ha lo, we can app ly the hydrodynamic ablation 
equations of lPittardI ()2007[ ) (pg. 245) to estimate the life- 
time of the model cloud. In this calculation, we use the 
total estimated cloud mass from Table [2] which includes 
all warm phase gas (WNM or WIM) and the lifetime is 
the time required for all of this warm phase gas to be 
heated and absorbed into the hot Halo. This lifetime es- 
timate will be on the high side since we are starting with 
the ablation process well underway with the warm enve- 
lope mass already in transition. It should be noted that 
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Fig. 9. — Neutral density characteristics of a 15% clumpy cloud core with uniform average density 5.5 X 10~^ cm~^, ionized 3D spherical 
cloud at Galactocentric distance 70 kpc exposed to standard Galactic and extragalactic ionizing radiation assuming the surrounding envelope 
is highly ionized and transparent. The cloud radius was adjusted to achieve the central column density of the cloud core observation-based 
model, Nhi = 4.8 X lO^^cm"^. Plots for each case are (a) LOS neutral column density; (b) density profile along a LOS through the 
centroid; and (c) average radial density profile along with that from our model. 



our averaging techniques expose a normally undetected 
envelope whose observed neutral mass is ~ 7 times that 
of the core, so a very large portion of the mass of this 
cloud is in the envelope. This suggests that ablation time 
calculations based on cloud masses inferred from column 
densities of conventional H i maps will be quite low since 
the unobserved envelope may contain the majority of un- 
ablated mass. 

Assuming a distance D ~ 70 kpc, the total cloud 
mass estimate from Table [2] is Afc = 4 x 10^ Mq. We 
assume a peculiar cloud velocity of 400 km s""'^ with 
respect to a non-rotating Halo (v^os — 350 km 
plus some moderate proper motion), WNM/WIM cloud 
sound speed Cc w lOkms"^, and Halo properties 
Thaio ~ lO^K, and Uhaio « lO-^cm-^ (iFang et al.l 
|2006[) on order of typically assumed values. The flow 
will be supersonic for these conditions, so we use M = 
{McCc)Hphaiov)i « 3 X lO^MeMyr"^ for a cloud life- 
time of tc ~ 140 Myr. The infall distance prior to de- 
struction is on order of vlos^c — 50 kpc, not quite sur- 
viving to the Galactic Disk from the assumed distance. 



However, mass loss rates could be considerably lower due 
to factors not accounted for in the above equation such 
as recycling of WNM/WIM gas through rap id cooling 
(iBegelman fc FabianI 119 90^ 'Kw aket al.l 120111 ). conduc- 
tion effects (iVieser fc He nslcr 2007a|) and magnetic fields 
lEsquivel et al."2006^. 

'Kwak et al. (2011) performed long term simulations of 
clouds with properties on the same order as ours except 
for an assumption of solar met allicity rather th an the 
0.1 solar estimate for this cloud (|Fox et al.ll201Cll ). Their 
model "D" has comparable initial mass and retains 70% 
of its mass after 240 Myr. Its velocity is less than half 
that of ours, but the simulations and the authors' analy- 
sis suggest that once supersonic, flow velocity has a large 
effect on cloud morphology, but little effect on the mass 
loss rate. These models suggest the cloud should have 
a lifetime as ^ yNM/^YIM on order of twice that calcu- 
lated with the iPittardI ()2007[ ) equation, long enough for 
it to potentially get within reach of the Galactic Disk. It 
should be pointed out that if cooling in the TML indeed 
slows ablation through recycling as mentioned above, the 
lower metallicity would reduce the lifetime by some un- 
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known amount. 

On the other hand, iHeitsch fc PutmanI (|2009() per- 
formed 3D simulations of clouds falling through the Halo 
assuming 0.1 solar abundances so cooling effects are ap- 
propriate to our cloud. They directly traced the smaller 
WNM/WIM clump fragments like those that compose 
our observed envelope. The cases illustrated are mostly 
subsonic at all times and clumps generally trail the core 
mass, but in one of the free-fall simulations shown, it 
briefly becomes transonic and develops some leading 
clumps suggesting that vortices are developi ng as t hey d o 
in other supersonic simulations like Kwak et al.l (j2011[ ). 
Their model "Wblal5b" has initial mass 3 x 10^ M© 
and loses 50% of its mass to the Halo in 50 Myr (their 
Figure 3) and extrapolating to full mass loss, lasting 
~ 100 Myr. To estimate our cloud's lifetime, we start 
with the integral core mass only (without the envelope) 
as this corresponds to the initial state of the simulations. 
We assume that the core will ablate similarly to the sim- 
ulated cloud independent of a the pre-existing envelope. 
Our neutral core mass at 70kpc is 2.5 x 10'^ M© (Ta- 
ble m and allowing for 24% ionization as estimated in 
Section [4.31 is ^ 3x the mass of the simulated cloud. 
The flow is MachO.7 in the simulation and so should ab- 
late at nearly the same rate. The extrapolated lifetime of 
our cloud would be on order of 300 Myr, more than twice 
the iPittardI ()2007f ) calculation, and easily sufficient to 
reach the Disk with some WNM/WIM mass remaining. 

In summary, the low-end estimate of cloud survival 
time as WNM/WIM prior to being subsumed by the Halo 
is 140 Myr, insufficient to reach t he Ga lactic Disk. This is 
based on the analysis ofl PittardI (|2007| ) which does not ac- 
count for cooling or other factors that have b een shown to 
suppr ess ablation. The numerical results of iKwak et al.l 
(|2011l ) which include cooling, suggest that the time to full 
destruction of warm gas in a similar cloud is much longer, 
by up to a factor of 2 or more although assumed solar 
metallicity puts this e stimate on the high s ide. H owever, 
rough extrapolation of'Heit sch fc PutmanI (|200 9') simula- 
tions, with abundances appr opriate to our cloud indicate 
that, at more than twice the IPittardI (|2007[ ) lifetime, it is 
possible or even likely that warm gas will reach the Disk. 

5. SUMMARY AND CONCLUSIONS 

We have obtained deep Hi 21 cm emission measure- 
ments mapping the vicinity of a relatively isolated, cir- 
cularly projected cloud in the northern MS and ap- 
plied a new spatial averaging technique to profile its 
mean spectrum vs. projected distance from the cloud 
center. Although near a larger cloud complex from 
which it appears to be separating, the cloud is well de- 
fined in angular extent and very well defined in veloc- 
ity. By exploiting the approximate azimuthal symme- 
try of the cloud and averaging along successively larger 
concentric circular paths, we obtained correspondingly 
increased sensitivity without compromising the essen- 
tial spatial resolution along the radial dimension. With 
this method, we detected and characterized core emis- 
sion having « 20 km s"^ WNM/WIM-like FWHM with 
velocity of « —340km s~^ and an envelope with broad 
diffuse emission having « 60 km FWHM and veloc- 
ity of ~ —360 km s^^ extending well beyond the appar- 
ent periphery of the cloud's emission. We then obtained 



robust estimates of the core and envelope spectral prop- 
erties by optimally matching the observed, averaged pro- 
file to a simple, spherical 3D core-plus-envelope emission 
model of the cloud projected on to the sky and averaged 
using the same method. 

The envelope probes the neutral component of a 
boundary layer between the infalling WNM/WIM cloud 
core and the ambient HIM of the Halo. Assuming that 
the broad line width reflects very warm transition gas in 
a conductive boundary layer, we estimate the detected 
extent of the neutral envelope is more than an order of 
magnitude beyond what should be detectable based on 
the analyses and simulations of cloud evaporation. We 
conclude that a Turbulent Mixing Layer (TML) best ex- 
plains the observed characteristics. Theoretical treat- 
ments and numerical simulations in the literature show 
that a TML has properties that can explain both the 
envelope line width and leading velocity where the en- 
velope is composed of small clumps of narrower-width 
WNM/WIM dispersed over a large velocity range amidst 
warmer gas in the TML. This presents a wide Gaussian- 
like spectrum when a large number or these small clumps 
are within the beam with some leading and some lag- 
ging the core component. The leading clumps are also 
explained by the TML where vortices develop near the 
infalling core when flow is supersonic as expected and 
demonstrated in various numerical simulations. 

Passing through the cloud and fortuitously near its 
center, is the NGC 7469 background source sight line 
which was characterized by FSWIO using optical and UV 
absorption data. Because O i and H i trace virtually iden- 
tical profiles weighted toward the center of a partially 
photoionized cloud, we used measured velocity disper- 
sions from FSWlO's O i and our H i core velocity compo- 
nents to establish a core temperature T — 8350 ± 350 K 
with a turbulent component at = 5.3km s~^, consistent 
with WNM/WIM. 

The role of photoionization on the cloud was investi- 
gated using a 3D photoionization equilibrium program. 
Assuming that the coUisionally ionized warm TML gas as 
well as the interspersed WNM/WIM clumps are highly 
ionized, we treated the envelope as transparent to Galac- 
tic and extragalactic ionizing radiation. We exposed 
a clumpy spherical cloud with radius 320 pc at a dis- 
tance of 70 kpc with uniform average Hydrogen density of 
5.5 X 10^"^ cm^"^ to combined Galactic and extragalactic 
ionizing flux based on published values at this distance. 
We determined that the centrally concentrated exponen- 
tial profile of the neutral component is approximated well 
by our simple uniform neutral density model of ^ 1/2 the 
radius. The core's actual extent as WNM/WIM is there- 
fore on order of twice that of our 3D neutral component 
model. The neutral fraction along the central LOS is 
24% and the ionization parameter is 9.7 x 10~^, quite 
different from FSWlO's estimates of 7% and > 5 x 10"^, 
respectively. We are unable to approach FSWlO's values 
with any realistic cloud at a realistic distance. We con- 
jecture that the discrepancy is at least partially due to 
FSWlO's inclusion of what we suggest are envelope veloc- 
ity components in abundance calculations. If produced 
in the clumpy, warmer and partially coUisional envelope 
environment we propose, these components would not be 
represented well by CLOUDY. 
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The cloud lifetime assuming simple hydrodynamic ab- 
lation driven by shear instability in supersonic flow is 
estimated at ~ 140 Myr, insufhcient for any WNM or 
WIM to reach the Disk and potentially fuel future star 
formation. However, a TML is expected to partially sup- 
pres s ablation through re-cooh ng of some of the heated 
gas ()Begelman fc Fabianlll99ClD. Simulati on of a compa- 
rable infalling cloud bv lKwak et all (|2011[ ) which includes 
cooling indicates a dramatically increased 2x lifetime, 
suggesting that a large portion of the total mass could 
reach the disk. T his result is s upported by extrapolating 
the simulations of iHeitsch fc Putmai] (j2009f ). 

Through our use of the spatial averaging technique, 
this work provides the first look (or more accurately, a 
peek) at the structure of the highly ionized and diffuse 
boundary layer between the warm neutral gas of the MS 
and the HIM of the Halo from the perspective of the 
low temperature side of the interface. FSWIO provides a 
rich set of data and analysis for the fortuitously aligned 
NGC 7469 sight line, creating a mini-laboratory to study 
MS-Halo interaction. We were able to incorporate some 
of their results to complement ours, but were unable to 
quantify characteristics of the barely probed TML. To 
fully exploit this alignment of complementary data, hy- 
drodynamic simulations that emulate the evolution of 
the larger environment of the cloud (i.e., including the 
neighboring cloud complex) would be needed that allow 
simulated observations on a similar cloud feature. Pho- 
toionization would have to be accounted for and its pos- 
sible impact on the hydrodynamic evolution would have 
to be addressed. 

Future work will include applying the technique to iso- 
lated circular or filamentary structures (with or without 



benefit of well-placed background sources) in other parts 
of the MS, both upstream and downstream to see if this 
apparent mixing layer is common and if its characteris- 
tics vary with time, indicated by the location along the 
MS. We also plan to obtain higher resolution maps on 
this cloud and others with the Arecibo Telescope to com- 
plement the lower resolution GET data and gain insight 
into the smaller scale structure. Further in the future, 
we plan to apply the technique to deep, very high resolu- 
tion observations of such clouds using a large array such 
as the Australian Square Kilometre Array Pathfinder 
(ASKAP). 
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